Strain engineering magnetic frustration in perovskite oxide thin films 
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Our first-principles results show that geometric frustration can be induced in thin films of multi- 
ferroic BiFeOs. We find that competing magnetic interactions occur in the so-called super-tetragonal 
phase of this material, which can be grown on strongly-compressive substrates. We show that the 
frustration level can be tuned by appropriately choosing the substrate; in fact, the three-dimensional 
spin order gets totally annihilated in a narrow range of epitaxial strains. We argue that the effects 
revealed here are not exclusive to BiFe03, and predict that they also occur in multiferroic BiCo03. 
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The competition between different interactions un- 
derlies some of the most fascinating phenomena in 
condensed-matter physics. When the competing orders 
are of similar strength, rich phase diagrams are likely to 
emerge, and the materials tend to be strongly responsive 
to external perturbations. Examples abound in the fam- 
ily of perovskite oxides, ranging from magnetoresistive 
manganites [l( to ferroelectric relaxors In addi- 

tion to the fascinating science they involve, competing 
interactions often lead to important functionalities. 

Of special interest are the cases in which the competi- 
tion relics on the topology of the crystal lattice, as in 
a triangular network of anti-ferromagnetically coupled 
spins, which is the classic example of geometric frustra- 
tion. Here we present first-principles results showing that 
strain engineering (i.e., taking advantage of the mismatch 
stress exerted by a substrate on a thin film) can be used 
to induce a novel kind of tunable geometric frustration in 
a spin system. Our work focuses on room-temperature 
multiferroic BiFeOs, and we show that the same effects 
also occur in other compounds like BiCo03. 

Tunable frustration in BiFeOs.- Multiferroic BiFeOs 
(BFO) is a perovskite oxide that presents ferroelectric 
(T c « 1100 K) and magnetic (T N « 760 K) orders at 
ambient conditions [4j. While the usual BFO phase is 
rhombohedral (R3c space group), it was recently dis- 
covered that BFO's atomic structure can be drastically 
modified by growing thin films on strongly compressive 
(OOl)-oriented substrates ||. The phases thus obtained 
are called super-tetragonal (ST), as their (pseudo-cubic) 
unit cell presents a very large aspect ratio c/a w 1.25 
(see Fig. 1). Such ST-BFO phases display a number of 
appealing features and are currently receiving a lot of at- 
tention; in particular, they undergo both structural and 
magnetic-ordering transitions slightly above room tem- 
perature (T r ) |6j,[7[, which might lead to improved func- 
tional properties. 

The magnetic order of ST-BFO remains an open prob- 
lem. Several first-principles works [H [§] predict the so- 
called C-type anti-ferromagnetic (C-AFM) spin arrange- 
ment, while most experimental studies suggest that the 
so-called G-AFM order dominates [l(j|. As shown 
in Fig. 1, the C-AFM and G-AFM orders are identi- 
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FIG. 1. (Color online.) Left: Elemental 5-atom cell of ST- 
BFO, indicating the in-plane epitaxial constraint a — b — 
a Bub . Right: Spin arrangements considered in this work. We 
sketch the 2x2x2 repetition of the elemental cell that we 
simulated, showing only the Fe atoms. In the FM case we 
indicate the exchange interactions J discussed in the text. 



cal within the ab plane (first-neighboring spins are anti- 
parallel), but differ along the out-of-plane z direction 
(first-neighboring spins are parallel in C-AFM and anti- 
parallel in G-AFM). It is generally accepted [6l.l8l-lllj that 
the in-plane exchange interaction between neighboring Fe 
atoms ( J a in Fig. 1) is anti-ferromagnetic and relatively 
strong; in contrast, the out-of-plane couplings (J c and 
J ac in Fig. 1) are believed to be small because of the 
very large separation between irons along the z direc- 
tion. In fact, MacDougall et al. [l(| have argued that 
the occurrence of C-AFM or G-AFM orders in specific 
films may be decided by factors that would be extrinsic 
to a perfect ST-BFO lattice, as the intrinsic out-of-plane 
couplings can be expected to be negligible. 

Wanting to shed light on these issues, we used first- 
principles methods [lj 



131 ] to investigate the magnetic 



order in ST-BFO as a function of the epitaxial strain ex- 
erted by a (OOl)-oriented square substrate. For this pur- 
pose, we considered a perfectly tetragonal atomic struc- 
ture (P4mm space group); we checked that, as regards 
the spin couplings, this structure is representative of the 
variety of lower-symmetry (monoclinic) phases [tj [13] 
that occur in the actual films. For each considered value 
of the substrate lattice parameter a sub [lit ], we stud- 




ated to a Fe 3+ -0 2 -Fe 3+ chain forming a 180° angle, for 
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FIG. 2. (Color online.) Panel (a): Energies of the different 
spin arrangements as a function of <z sub . Panel (b): Exchange 
constants defined in Fig. 1 as a function of a Buh . J ac is mul- 
tiplied by 4, to better visualize the point at which Eg ~ Ec- 



ied several spin arrangements (FM, A-AFM, G-AFM, 
and G-AFM, sketched in Fig. 1); for each arrangement, 
we relaxed the atomic structure subject to the epitaxial 
constraint. As shown in Fig. 2(a), we find that the C- 
AFM arrangement is the lowest-energy solution for the 
most stable films, which correspond to a snh ~ 3.71 A. 
The G-AFM order gets stabilized as a sub increases, the 
transition between C-AFM and G-AFM occurring at 
a sub » 3.79 A. 

To gain more insight, let us consider the Heisenberg 
spin Hamiltonian E — E = 1/2 N y^ytj JijSi • Sj, where 
N is the total number of Fe atoms, needed to capture 
such a crossover. Let us assume classical spins and take 
| Si | = 1 for simplicity. Then, if we restrict ourselves to 
couplings between first-nearest neighbors (i.e., to J a and 
J c in Fig. 1, as discussed by other authors (a. Hoj). the 
energies per Fe atom of the C-AFM and G-AFM orders 
are E c = E Q - 2J a — J c and E G = E — 2J a + J c , re- 
spectively. The crossover point (Eg = Ec) would thus 
correspond to J c — 0, i.e., to the value of a suh at which J c 
changes sign. However, this simple model is not satisfac- 
tory, for two main reasons: (1) It is difficult to imagine 
an exchange mechanism that may render a ferromagnetic 
(FM) coupling J c < 0. Indeed, the J c coupling is associ- 



which the well-known Goodenough-Kanamori rules 16 1 
predict an AFM super-exchange interaction; such an ex- 
pectation should hold even in a case like this one, where 
the two Fe-0 bonds in the super-exchange path have dif- 
ferent lengths. (2) For J c < 0, such a simple model pre- 
dicts Ec < Eg < Ep < Ex- However, our results show 
that the FM order is the least favorable one through- 
out the considered a suh range. Indeed, our calculations 
predict that, even though there is a crossover between 
the C-AFM and G-AFM orders, where the magnetic in- 
teraction between ab layers changes sign, the analogous 
crossover between FM and A-AFM does not take place! 

These difficulties are resolved by extending the model 
in the simplest possible way, i.e., by including the J ac 
interaction defined in Fig. 1. The energies per Fe atom 
of our four magnetic orders are then given by: 



E-p — Eq + 2J a + J c + 4 J ac , 
Ea = E + 2J a — J c — AJ ac , 
Ec=E - 2J a + J C - 4J ac , 
E G =E - 2J a -J c + AJ ac . 
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We can thus compute the J's from the data in Fig. 2(a); 
the results are shown in Fig. 2(b). This extended model 
is able to reproduce exactly the first-principles energies 
of our four magnetic structures in the whole a suh range. 
In addition, the obtained values of J c and J ac are always 
positive, which corresponds to AFM interactions; as re- 
gards </ c , this is compatible with the expectations from 
the Goodenough-Kanamori analysis. 

It may seem surprising that this model can capture the 
C-AFM ground state, given that the computed out-of- 
plane interactions, J c and J ac , are both AFM in nature. 
To understand this, let us consider first what determines 
the relative stability of the FM and A-AFM orders, for 
which we have E-p — Ea = 2J C + 8J ac - Here, in both 
cases the ab planes are ferromagnetically ordered, and 
both J c and J ac favor the A-AFM solution over the FM 
one. On the other hand, the energy gap between C-AFM 
and G-AFM is E c ~ E G = 2J C - 8J ac . In this case, 
the ab planes are anti-ferromagnetically ordered, and the 
out-of-plane interactions compete: A positive J c favors 
the G-AFM order via the AFM coupling between Fe ions 
that are first neighbors out-of-plane. In contrast, a pos- 
itive Jac favors the C-AFM order via the AFM coupling 
between Fe ions that are second neighbors out-of-plane. 
When J ac is large enough - i.e., when 4J QC > J c , where 
the factor of 4 comes from the ratio between first- and 
second-nearest neighbors out-of-plane-, the C-AFM or- 
der prevails. Hence, according to this model, the pre- 
dicted stabilization of the C-AFM phase relies on both 
the AFM interaction J ac and the presence of a robust 
AFM order within the ab planes. 

The obtained a sub -dependence of the exchange con- 
stants J a and J c seems rather natural. As we compress 




FIG. 3. (Color online.) Results from Monte Carlo simulations. Panels (a) and (b): T-dependence of the order parameters 
quantifying the degree of C-AFM (Ac) and G-AFM (^4g) order, respectively. We show the results for all the substrates 
investigated, with a sub ranging from 3.65 A to 3.84 A. A s = iV" 1 Sj Z exp (iq a Rj), where N is the number of cells in 
the simulation box and Rj is the lattice vector corresponding to spin Sj; q s defines the s-like spin arrangement, with qc = 
27r/a sub (l/2, 1/2, 0) and q G = 27r(l/2a sub , l/2a sub , l/2c). We onlyneed to consider the z component of the spins (S jz ) because 
of the small symmetry-breaking included in our Hamiltonians 17]. A perfect s-like order corresponds to having A s — 1. 
Panel (c): Ac and Ac obtained at 50 K and as a function of a sub . Panel (d): results for a sub = 3.782 A. Ai aye r quantifies the 
AFM order within the ab layers. Ai ayor = N{^ yet • Sj z exp (iqzoRj), where the primed sum runs over the spins in the first ab 
layer, which is representative of the rest; Ni ayet is the number of cells in a layer and q2D = 2-7r/a sub (l/2, 1/2, 0). 



in-plane, the Fe spins coupled by J a get closer and their 
interaction becomes stronger; in contrast, the distance 
between irons coupled by J r grows (results in |13|) and 
the interaction weakens [18j]. However, it is not clear 
what to expect for J ac . In this case, any coupling mech- 
anism that one can imagine will have both in-plane and 
out-of-plane components; hence, it is probably not sur- 
prising to find that J ac varies weakly with a sub . 

Hence, our calculations and model analysis reveal a ro- 
bust mechanism leading to a crossover between G-AFM 
and C-AFM orders as a sub varies. Such a transition is 
the result of the competition between two magnetic in- 
teractions that become comparable in a certain range of 
epitaxial strains. In fact, our ST-BFO films can be con- 
sidered a case of frustrated spin system, where the mag- 
nitude of the frustration can be tuned by appropriately 
choosing the substrates on which the films are grown. 

Phase diagram.- We solved our Hamiltonians by 
performing Monte Carlo simulations in a periodically- 
repeated box of 20x20x20 spins [H, HJ. Figures 3(a) 
and 3(b) show, respectively, our results for the T- 
dependence of the parameters monitoring the C-AFM 
(A c ) and G-AFM (A G ) orders; Fig. 3(c) shows the re- 
sults for T = 50 K, where the evolution of the ground 
state with a sub is easily appreciated (l7| . 

Our simulations render a paramagnetic (PM) phase at 
high temperatures. Then, for a sub < 3.778 A we observe 
a transition to a C-AFM phase as T decreases; in con- 
trast, for a sub > 3.792 A the low-T phase presents the 
G-AFM spin order. In the narrow intermediate region 
3.778 A < a sub < 3.792 A, some sort of T-driven tran- 



sition occurs as evidenced by the non-zero values of Ac 
and Ac', however, no clear-cut three-dimensional order 
emerges. 

Figure 3(d) shows results for a sub = 3.782 A, which is 
representative of the intermediate region. In this case, 
Aq and Ac take values that look rather arbitrary. Nev- 
ertheless, if we consider the AFM order within the ab 
planes, as quantified by Ai aycr defined in the figure cap- 
tion, we recover a well-behaved transition. The picture 
that emerges is thus clear: For intermediate values of 
a sub ^ ne S y S t em becomes two-dimensional (2D). The out- 
of-plane correlations are very weak, and long-range order 
along z essentially disappears. 

Figure 4 shows the phase diagram that emerges from 
our results. Note that the width of the 2D region depends 
on the specific box size and sampling method used in 
our simulations. Indeed, for small but non-zero values of 
Jc — 4 J ac | ) the ground state of the system is well denned 
in the thermodynamic limit. However, the specifics of 
our simulations will determine whether the AFM planes 
can order correctly along the z direction or, instead, get 
stuck in a meta-stable configuration displaying disorder 
and/or phase co-existence. [As appreciated in Fig. 3(c), 
we tend to find incomplete G-like (resp. C-like) order on 
the right (left) side of the intermediate region, which is 
a result of having J c — AJ ac > (resp. J c — 4J ac < 0).] 
Hence, the 2D region in Fig. 4 has to be taken as evidence 
for what probably is a line separating the C-AFM and 
G-AFM phases in the ideal case. 

In order to understand the evolution of Tn with a sub , 
it is useful to resort to a mean-field (MF) analysis of our 
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FIG. 4. (Color online.) Phase diagram deduced from our 
Monte Carlo simulation. The dashed lines show the transition 
lines obtained at the mean-field level, which are rescaled to 
fit with the Monte Carlo results (see text). 



Hamiltonians. Within this approximation, the transition 
temperature is proportional to the mean field experienced 
by the spins in the lattice; such a field is essentially given 
by 4J Q - 2 J c + 8J ac in the C-AFM case, and by 4J Q + 
2 J c - 8J ac for the G-AFM solution. In Fig. 4 we plot the 
a sub -dependence of these MFs; note that we multiply the 
results by a factor that is constant through the whole a sub 
range and common to both spin arrangements, so that 
the MF lines fit in the scale of the Monte Carlos results. 
As we can see, the general dependence of Tn with a sub 
is captured at the MF level. In particular, for decreasing 
a sub values, T N of the PM->C-AFM transition grows as 
a consequence of the increasing in-plane coupling J a . On 
the other hand, if a sub becomes larger, Tn of the PM^G- 
AFM transition grows driven by the rapidly increasing 
J c interaction. Finally, as we approach the 2D region, 
the Tn values obtained from Monte Carlo decrease more 
strongly than we would predict at the MF level. This 
is most likely a consequence of the competition at play 
in the Monte Carlo simulations, as such an effect is not 
captured by our simple MF analysis. 

Further discussion.- Previous works have shown that 
the approach adopted here renders reliable Neel temper- 
atures for ST-BFO Indeed, for BFO films grown 
on LaAlC>3 (a suh = 3.79 A), measured Tn values range 
from 324 K [lOj to 360 K 16 1 , which is in reasonable agree- 
ment with our results [2l|, |22j . 

Regarding the specific magnetic order, as far as we 
know only films grown on LaAlOa have been charac- 
terized experimentally. There seems to be consensus 
about the fact that the films adopt a G-AFM order at 
T r (EH, EH- Additionally, MacDougall et ai. have 



reported a coexistence of G-AFM and C-AFM orders 
at temperatures below 260 K. Such a situation seems 
compatible with the 2D region that we find precisely 
for a sub w 3.79 A. Indeed, for J c - 4J ac > we ob- 
tained a dominant, but not complete, G-AFM order [see 
Fig. 3(c)], which clearly resembles the experimental find- 
ings of Ref . 

At and around the 2D region we have J c — AJ ac w 0, 
which results in an effective magnetic decoupling of the 
ab planes. In such cases, it is conceivable that impuri- 
ties, defects, and other extrinsic factors may influence 
the nature of the magnetic ground state or the relative 
populations of co-existing phases. However, for smaller 
a sub vames we obtain a sizable J c — 4 J ac < 0. Hence, we 
predict that C-AFM order will occur in ST-BFO grown 
on substrates more compressive than LaAlOa. 

Note that this picture differs significantly from that 
of MacDougall et al. [l(J, who argued that all the out- 
of-plane couplings are essentially negligible (i.e., J c « 
and J ac « 0) in ST-BFO grown on LaA10 3 . Their inter- 
pretation implies that the magnetic order in ST-BFO 
should be 2D-like for any substrate more compressive 
than LaA103, and one would not expect to observe any 
robust C-AFM phase. This is clearly at variance with 
our results. 

How general is the tunable geometric frustration pre- 
dicted in this work? Let us emphasize that the require- 
ments to observe these effects do not seem exotic at all: 
In essence, we need to (1) have AFM interactions between 
neighboring magnetic cations, and (2) obtain an expan- 
sion of the out-of-plane lattice parameter as we compress 
in-plane. Condition (1) is satisfied by many perovskite 
oxides that are AFM insulators. Condition (2) is the 
expected behavior for all perovskites; further, there is a 
growing number of compounds that are known to dis- 
play large c/a ratios when grown under appropriate con- 
ditions. 

To test this presumed generality, we investigated 
BiCoOa (BCO) under epitaxial strain. BCO is an insu- 
lating perovskite whose ground state has the ideal ST 
structure {PAmm space group, with c/a = 1.27 and 
a =3.73 A) and a C-AFM spin order [23j. As we im- 
posed a tensile epitaxial strain, we observed a transition 
from C-AFM to G-AFM at a sub « 3.84 A. As in the case 
of ST-BFO, the analogous crossover between the FM and 
A- AFM orders is absent, and the transition point is de- 
fined by the condition J c = 4J ac , with all the computed 
exchange couplings being AFM in nature. Hence, BCO 
presents exactly the same magnetic frustration effects 
that we have discussed for ST-BFO. Our BCO results 
are summarized in [l3j |. including the prediction that G- 
AFM order can be obtained in BCO films at T r . Hence, 
BCO, as well as the BFO-BCO solid solutions studied in 



Ref. |!Jj, may offer interesting alternatives for the experi- 
mental investigation of these effects. 
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